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ABSTRACT 

Due to the short settling times of metals in DA white dwarf atmospheres, any white dwarfs with photospheric 
metals must be actively accreting. It is therefore natural to expect that the metals may not be deposited uni- 
formly on the surface of the star. We present calculations showing how the temperature variations associated 
with white dwarf pulsations lead to an observable diagnostic of the surface metal distribution, and we show 
what constraints current data sets are able to provide. We also investigate the effect that time-variable accre- 
tion has on the metal abundances of different species, and we show how this can lead to constraints on the 
gravitational settling times. 

Subject headings: accretion, accretion disks — diffusion — convection — stars: oscillations — stars: variables: 
other — white dwarfs 



1. ASTROPHYSICAL CONTEXT 

There are two main classes of white dwarf stars: those with 
hydrogen-rich atmospheres (spectral type DA) and those with 
helium-rich atmospheres (non-DA spectral types). The rea- 
son for this is that the high surface gravities of white dwarfs 
lead to efficient gravitational settling, with the lightest ele- 
ments rising to the surface. In addition, some DA white 
dwarfs have spectra showing metal l ines of elements such as 
Ca and Mg (Zucker man et al.l 120031) . and these stars are re- 
ferred to as DA Zs; about 20% of a ll DAs fall into this cate- 
gory. Recently, iDufour et al.l (120071) have announced a new 
class of white dwarf with carbon-dominated atmospheres, the 
"hot DQ" stars, several exampl es of which have b een found 
in the Sloan Digital Sky Survey dLiebert et al.ll2003l) . 

The presence of metals in the DAZs is intriguing since 
the settling time scale for the metals may be many or- 
ders of magnitude shorter than the evolutionary age of 
these objects. Indeed, for DAZs with T e g ~ 12,000 K, 
the settling time scale can be on the order of days or 
weeks, meaning that these object s are experiencing on- 
going accretion dKoester & Wilkenl 120061) . This ongoing 
accretion is consistent with the fact that nearly a dozen 
of t hese objects have detected dust disks ( " debris disks," 
see iTokunaga. Becklin. & Zuckerrrlanl 119901: iBecklin et all 
20051: iKilic et al 1 120051: iReach et all 120051: iKilic et al.l | 2006t 
von Hippel et al.l 120071: lJura . Farihi, & Zucker man! 120071: 
Juraet al.1 120071 iFarihi. Zuckerman. & Becklinl 120081) "and 
these disks are assumed to be the sources of the metal lines 
seen in these white dwarf atmospheres. The best studied 
object of this DAZ class with an observed disk, G29-38, is 
also a multi-periodic variable white dwarf (DAV), pulsating 
in non-radial g-modes with periods of a few hundred to one 
thousand seconds. 

The technique of asteroseismology uses the observed pul- 
sation modes of a star to infer and constrain the interior 
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structure of the star, thus obtaining information on the star's 
structure as a function of radius (Bradley & Winget 1994; 
[ Kawaler & Bradlevl[T994l: iMetcalfe. Salads, & Wingetll2002 : 
Montg omery. Metcalfe. & Wingetll2003l) . In contrast to this, 
our approach in this paper is to use the different angular de- 
pendence of the pulsations to constrain the accretion process. 
Since the accretion is most likely occurring through a disk, 
it is natural to suppose that the metals may not be uniformly 
distributed across the star's surface. We present calculations 
showing how we can place constraints on the non-uniformity 
of the accretion process. 

2. GRAVITATIONAL SETTLING AND HORIZONTAL DIFFUSION 

For all our calculations, we use the DAZ G29-38 as our 
template, since, given its brightness, it has the greatest po- 
tential for successful measurement of a surface inhomogene- 
ity. In addition, we have archival data on this star appropriate 
to this application. In this section, we will therefore attempt 
to quantify the importance of gravitational settling and hori- 
zontal diffusion assuming a model with parameters similar to 
those of G29-38. 
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FIG. 1 . — Possible surface metal distributions for accretion centered on the 
poles (dashed curve) or the equator (solid curve) as a function of the polar 
angle (co-latitude). 
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FIG. 2. — The diagnostic R as a function of inclination angle The left panel shows the results for the polar distribution of metals shown in FigurefT]and the 
right panel shows the results for the equatorial distribution. The different curves are labeled by the 1 1 m) values of the relevant pulsation modes. 



iBergeron et all d2004l) find T eS = 1 1,820 K and l og,? = 8.14 
for this star, while lKoester. Provencal. & Shipmanl(ll997l) find 
r e f t = 1 1,600 K and logg = 8.05. Interpolating in the tables of 
Koester & Wilken for Ca yields a settling time of ~ 1 3 days 
for the first set of parameters 5 and a settling time of ~ 23 days 
for the second. We note that while unseen helium in the atmo- 
sphere could lengthen th ese settling times considerably (e.g., 
iGarcfa-Berro et al. 2007), its presence is inconsistent with the 
pulsation results for this star: such an amount would imply 
a much deeper surface conve ction zone, in conflict with that 
found bv lMontgomervl d2005l) . 

Since these stars should have surface convection zones, the 
dominant form of horizontal transport of Ca will be due to 
the turbulent viscosity. We can estimate the size of this dif- 
fusion coefficient as D ~ Vc4, where Vc is a typical con- 
vective velocity and /;, is the assumed "mixing length" for 
convection. From our w hite dwarf evolution code (e.g., see 
Montg omery et al] [l999) we find for both sets of stellar pa- 
rameters that D w 1.5 • 10 I0 cm 2 /sec. 

For these simple exploratory calculations we assume az- 
imuthal symmetry for both the accretion and the surface metal 
distribution, i.e., Z = Z(8, t) and S = S(8, t), where Z(8, t) is the 
metal abundance, S(8,t) is the source function of metals ac- 
creting onto the white dwarf, 8 is the "co-latitude" of a point 
on the star's surface, and t is time. 

Since convection will uniformly mix material vertically, we 
can treat the convective region as a single "zone" and write an 
equation for the time evolution of Z as a function of 8 and f : 



dZ(8,t) 
dt 



-~fZ(8,t)+DVlZ(8,t) + S(8,t), 



(1) 



where 7=1/ (settling time) is the settling rate, V{; is the hor- 
izontal part of the Laplacian operator, and the other vari- 
ables are as defined above. An estimate of the relative 
importance of sinking to spreading is r\ = jR^/D, where 
R* is the radius of the white dwarf: r\ ~ 50 using the 
Bergeron et alj (120041) values wherea s 77 ~ 30 for those of 
Koester. Provencal. & Shipmanl (119971) . When 77 ^> 1 the met- 
als will sink before they have a chance to diffuse horizontally, 
while for 77 <C 1 the metals will have a chance to mix thor- 



Due to a mis- rea ding of Table 2 in Koester & Wilken 12006). 
von Hippel & Thompson (2007) calculated a settling time of ~ 7 days. 



oughly horizontally before sinking, producing a nearly uni- 
form surface distribution. 

In Figure Q] we show the metal distributions which arise 
from solutions of equationQ] The dashed curve is the equilib- 
rium distribution which results from constant accretion cen- 
tered at the poles and the solid curve is that which results from 
constant accretion centered on the equator. 

3. THE DIAGNOSTIC 

The flux variations observed in pulsating white dwarfs are 
due almost entirely to temperature changes o n the surface 
of the stars dRobinson. Kepler. & Nathenll982l) . These same 
temperature changes will also affect the equivalent widths 
(EWs) of any spectral lines, and in particular the EWs of metal 
lines. These metals may not be uniformly distributed across 
the star's surface and since the temperature variations are also 
non-uniform, we hope to be able to constrain the surface metal 
distribution. 

The relevant diagnostic we have developed, denoted by R, 
is the ratio of the fractional EW amplitude to the fractional 
flux amplitude: 



R 



S(EW)/(EW) 
6 F x I F x 



IFlux 



(2) 



where A EW is the fractional amplitude of EW variations, and 
A F i ux is the amplitude of the fractional flux variations observed 
in the given passband X. 

To compute R, we assume a particular surface tempera- 
ture perturbation of the form ST /T oc Yi m (8, </>). We then use 
model atmospheres to turn this into EW and flux variations on 
the surface of the star, taking into account the fact that the EW 
of the Ca lines will be directly proportional to the local abun- 
dance of Ca (e.g., the curves in Figure [TJ. For the passband 
X we assume a wavelength response appropriate to the Argos 
CCD with a BG 40 filter on the 2.1m teles cope at McDon- 
ald Observatory (Nather & Mukadam 2004); the wavelength 
range is taken to be 3000 A to 7000 A, with a peak response at 
5400 A. Finally, we integrate the result across the visible sur- 
face of the star. Limb darkening is automatically taken into 
account by this procedure. 

We note for a uniform distribution of metals that R = 2.1\ 
for I = 1 and R = 2.59 for 1 = 1. Thus, even by measuring R 
for a single mode it may be possible to tell whether metals 



Metal Inhomogeneity and Settling in Pulsating WDs 



3 



6 

< 

< 4 

II 

OS 



614.2 s 817.7 s 653.1 s 776.0 s 283.4 s 429.6 s 



uniform 



2 4 8 

mode number (dummy label) 

FIG. 3. — Values of R for modes in the DAV G29-38 as obtained from 4 
hours of Keck data in 1999. The values for the first 6 highest amplitude modes 
are shown, with the period in seconds indicated above each point. Only the 
first two modes have a signal-to-noise ratio large enough to be useful. 

are uniformly distributed on the star's surface. To obtain fur- 
ther constraints, we need R determinations from other modes 
and/or additional information such as the inclination angle of 
the star. 

In Figure [2] we show the R diagnostic as a function of in- 
clination angle 9\. The different curves are labeled according 
to their I and m values as \£m). The left-hand plot is for the 
polar distribution of metals given in Figure [T] and the right- 
hand plot is for the equatorial distribution. The thin horizontal 
boxed region shows the value of R expected if the metal dis- 
tribution is uniform; we see that a non-uniform distribution is 
very unlikely to produce a value of R in this range. 

4. COMPARISON WITH OBSERVATIONS 

In 1996 Clemens & van Kerkwijk obtained 
over 4 hours of time-resolved spectroscopy o f the 
DAV G29-38 dclemens. van Kerkwijk. & Wul 120001 
Ivan Kerkwijk. Clemens. & Wul l2000t IClemens etail 1 19991) . 
While not originally intended for this purpose, we can use this 
as an example data set for the technique proposed in the pre- 
vious section. We tak e the amplitudes for the EW vari ations 
from the analysis of Ivon Hippel & Thompson! (120071) and 
we use the amplitudes of the broadband (5200-5500 A) flux 
variati ons as determined by Ivan Kerkwijk. Clemens. & Wul 
(2000) 6 . We calculate R = Aew/Ahux for each of the modes, 
taking into account the errors on all quantities. 

We show the results of this procedure in Figure [3] Of the 
6 modes we identified, only the first two have R values with 
small enough error bars to provide any constraint on the Ca 
distribution. The 614 s mode provides the most convincing 
evidence for a non-uniform Ca distribution, since its R value is 
a full 3(7 above the range produced by a uniform distribution. 
The 818 s mode may also provide some evidence, although it 
is only about 1.3er above the value expected from a uniform 
distribution. 

With only one statistically significant value of R, we 
cannot hope to infer anything further about the Ca dis- 
tribution. Fortunately, we do have additional informa- 

6 Even though this wavelength range is much narrower than that used in the 
previous section for Figure [2] the central wavelengths of the two passbands 
are nearly the same and the results obtained are virtually indistinguishable 
from one another. 
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FIG. 4. — Comparison of the observed R value for the 614 s mode in G29- 
38 (shaded region) to the value expected for a polar, an equatorial, and a 
uniform distribution of Ca, as a function of inclination angle 8[ . The data are 
only consistent with the equatorial case. 



tion. IClemens. van Kerkwijk. & Wul (120001) used time- 
resolved sp ectroscopy to deter mine that this mode has £ = 1 . 
In addition, Montgomery (2005j), by modeling the nonlinear 
pulse shape of this mode, was able to further specify that it is 
an £ = 1, m = 1 mode. With just this one additional constraint, 
we are somewhat surprisingly able to constrain our models of 
the Ca distribution. FigureUshows the R values derived using 
an equatorial, a polar, and a uniform distribution, for a variety 
of inclination angles. The 614 s measurement is inconsistent 
with a polar distribution of metals, and is most consistent with 
an equatorial distribution. 

5. TIME- VARIABLE ACCRETION AND SETTLING 

In the case of G29-38, accretion is thought to be from a de- 
bris disk having roughly asteroidal composition, so that sev- 
eral elements are being accreted simultaneously. In particular, 
if we look at Ca and Mg, these two elements have different 
settling rates, with Mg settling about 50% more slowly than 
Ca. 

Accr etion is not expected to be a c ompletely steady pro- 
cess. Ivon Hippel & T hompson! d2007l) found the EW of Ca 
lines in G29-38 to vary on a time scale of weeks to years, 
with additional evidence for variations on shorter time scales 
(von Hippel, private communi cation). On the other hand, 
iDebes & Lopez -Morales! ( 120081) found the observed EW vari- 
ations of G29-38 to be consistent within the measurment er- 
rors. For the present we treat this as an open question and ask 
what effect time variable accretion would have on the mea- 
sured atmospheric abundances of different chemical species, 
and how this can be used as a probe of the settling process in 
these stars. 

For this problem it is sufficient to consider the spherically 
averaged version of equation [1] 



dZjjt) 
dt 



-jiZM+XiSd), 



(3) 



where Z, is the abundance of element i, 7* is its settling rate, 
and X, (assumed to be constant) is the relative fraction of ele- 
ment i in the accreted material. 

In Figure [5] we illustrate the solution of equation [3] for a 
particular set of parameters and a given accretion rate. The 
accretion rate assumed is shown in the top panel of Figure [5] 
For the two elements we have assumed that 71 = 1/13 d -1 and 
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FIG. 5. — Bottom panel: the fractional relative abundance of element 1 
(71 = 1/13 d~', solid line) and element 2 (72 = 1/26 d _1 , dashed line) as 
given by the solution of equation [5] Top panel: the time-variable accretion 
rate (arbitrary units) assumed for these calculations. 

71 = 1/26 d -1 . In the bottom panel we plot the abundance of 
each element divided by its mean value, i.e., its "fractional" 
abundance: the solid curve is that of element 1 and the dashed 
curve is that of element 2. By using relative/fractional abun- 
dances, we have divided out the effect of the X, parameter. 

The first feature to note is that the abundances of the two el- 
ements are highly correlated, i.e., when one goes up or down 



so does the other. In addition, we see that the amplitude of 
the excursions of element 1 are larger than those of element 2. 
In fact, taking the standard deviations of each of the curves 
we find that that of element 1 is about 1.8 times that of ele- 
ment 2. In other words, the element with the shorter settling 
time experiences larger fractional variations, and these vari- 
ations may be used for a rough estimate of the ratio of the 
settling times. Thus, periodic observations of a white dwarf 
with more than one metal line may be able to place constraints 
on the relative settling times of different chemical species in 
white dwarf atmospheres. 

6. CONCLUSIONS 

We have shown how the temperature variations due to stel- 
lar pulsation can be used to constrain the metal distribu- 
tion on the surface of a white dwarf, and we have shown 
that data currently in hand for the star G29-38 suggests that 
the metal distribution, and therefore the accretion, may be 
equatorial. If further studies support an equatorial distribu- 
tion for Ca in G29-38, this argues against magnetic accre- 
tion onto spots near the poles. Non-magnetic equatorial ac- 
cretion would further imply that the inner edge of the disk 
is much thinner than the white dwarf's radius. A physi- 
cally thin disk is consistent with observations to date (e.g., 
Uura. Farihi. & Zuckermanll2007t IvonHippel et al.ll2007l) . 

We have also shown how the observed variations in atmo- 
spheric abundance of two chemical species can be used to 
place constraints on their settling rates. The element with the 
faster settling rate will show the larger fluctuations, and, at 
least within the context of this simplified model, the ratio of 
the magnitude of these excursions from the mean will be ap- 
proximately proportional to the settling rate. 

Thus, the time variability of these systems, both in terms of 
the pulsations of the white dwarf and in terms of the variable 
accretion rate, allows us to probe the physics of accretion and 
gravitational settling in these objects. 
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